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We report time-of-flight neutron scattering measurements of the magnetic spectrum of Tb 3+ in 
Tb2Ti207. The data, which extend up to 120 meV and have calibrated intensity, enable us to 
consolidate and extend previous studies of the single-ion crystal field spectrum. We successfully 
refine a model for the crystal field potential in Tb2Ti2C>7 without relying on data from other rare 
earth titanate pyrochlores, and we confirm that the ground state is a non-Kramers doublet with 
predominantly [ ± 4 ) components. We compare the model critically with earlier models. 

PACS numbers: 75.10.Kt, 75.40.Gb, 71.70.Ch, 78.70.Nx 


I. INTRODUCTION 

Among the many magnetically frustrated pyrochlore 
oxides, Tb 2 Ti 2 07 (TTO) stands out because of its in¬ 
triguing low temperature state which is thought to be a 
type of spin liquid. 1 TTO shows no sign of any conven¬ 
tional symmetry-breaking transition (magnetic or struc¬ 
tural) down to temperatures as low as 50 mK (Refs, d 
and y) despite an antiferromagnetic Curie-Weiss tem¬ 
perature of —19 K (Ref. H) and predictions of magnetic 
order at 1 -2 K (Refs. d and 6). There are, however, strong 
short-range antiferromagnetic correlations in TTO at low 
temperatures^ - — 

The spin liquid state in TTO is not fully understood, 
but a key factor is the low energy part of the crystal- 
field-split / electron manifold of Tb 3+ , which comprises 
two doublets separated by about 1.5 meV (Refs. |2j and 
@). This splitting is comparable with the exchange and 
dipolar coupling strengths in TTO. Therefore, although 
the single-ion ground state is Ising-like, transverse fluc¬ 
tuations of the ground state moment appear once inter¬ 
actions are taken into account, and the cooperative na¬ 
ture of these fluctuations could lead to a quantum spin 
ice state i 11 ! 12 Moreover, the two doublets are also con¬ 
nected by quadrupolar interactions, allowing coupling to 
phonons. Indeed, the recent observation of hybridization 
between an acoustic phonon and a crystal field excitation 
in TTO demonstrates the importance of magnetoelastic 
interactions in this system^^ So the conditions exist in 
TTO for enhanced quantum fluctuations and suppressed 
magnetic ordering via both magnetic and magnetoelastic 
interactions, and an accurate determination of the crystal 
field states is required for quantitative modelling. 

Several models for the crystal field potential in TTO 
have been reported based on analyses of neutron and 
optical spectra combined with thermal and magnetic 
data lii 14 - 2 - Unfortunately, there are significant discrep¬ 
ancies between the published sets of crystal field param¬ 
eters. These are partly caused by the fact that some of 
the models employ a truncated basis containing only the 


13 states of the ground state 7 Fg manifold, whereas oth¬ 
ers include all the states in the / 8 configuration of Tb 3+ . 
The main difficulty, however, has been that up to now 
only four transitions within the J = 6 level have been 
observed unambiguously. Information from other heavy 
rare-earth titanates has generally been used to augment 
the data on TTO so that the six free parameters needed 
to describe the single-ion crystal-field Hamiltonian can 
be determined independently. 

In a very recent neutron scattering study^ a split¬ 
ting of one of the peaks was resolved and assumed to 
originate from two distinct crystal field transitions. This 
led to a significantly different crystal field potential com¬ 
pared with previous models for TTO. For instance, the 
wave functions of the ground and first excited doublets 
are interchanged compared with models informed by data 
from other rare-earth titanates. The correctness of this 
assumption has recently been questioned, and the split¬ 
ting attributed instead to coupling to a phonon^. 

The purpose of this paper is to resolve these discrepan¬ 
cies as far as possible. We report neutron scattering mea¬ 
surements of TTO which extend to higher energies than 
before. We have observed a crystal field transition that 
was not previously detected, and we have determined the 
transition intensities on an absolute scale. By fitting the 
data from TTO alone, using a model that includes in¬ 
termediate coupling basis states and partial J-mixing, 
we find a crystal field potential that provides a good de¬ 
scription of the experimental data. The model described 
in Ref. HO cannot be reconciled with our data. 


II. EXPERIMENTAL DETAILS 

A powder sample of mass 16 g was prepared by stan¬ 
dard solid state synthesis. The sample was found to be 
single phase and of high quality to within the precision 
of laboratory X-ray diffraction. Neutron inelastic scatter¬ 
ing measurements were performed on the MERLIN time- 
of-flight spectrometer at the ISIS Facility^ The sample 
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Figure 1. (Color online) Neutron scattering spectrum of poly¬ 
crystalline Tb 2 Ti 207 at T — 7K, measured on MERLIN with 
Ei = 150 meV. The colors represent the value of E x S(Q,E) 
in units of mb sr -1 f.u. -1 


was contained in an aluminium foil packet in the form 
of an annulus of diameter 40 mm and height 45 mm and 
sealed in an aluminium can containing helium exchange 
gas. The can was cooled by a closed-cycle refrigerator. 
Spectra were recorded for approximately 4 hours each 
with neutrons of incident energy Ei = 65 and 150 meV 
at temperatures of T = 7 and 90 K. The raw data were 
corrected for detector efficiency, sample attenuation and 
time-independent background following standard proce¬ 
dures. Vanadium spectra recorded at the same two in¬ 
cident energies were used to determine the energy reso¬ 
lution and to convert the intensities into units of cross 
section, mbsr -1 meV -1 f.u. -1 , where f.u. stands for the 
formula unit of Tb 2 Ti 2 07 . 

III. RESULTS 
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Figure [T] provides an overview of the data recorded at 
T = 7K with Ei = 150 meV. The spectrum is presented 
in the form of a color map of E x S(Q, E), where S(Q, E) 
is the scattering intensity as a function of energy E and 
the magnitude of the scattering vector Q = |Q|. Multi¬ 
plication by energy suppresses the strong elastic and low 
energy inelastic scattering and makes the weaker signals 
at higher energies more visible. 

Previous neutron measurements of the spectrum of 
TTO revealed magnetic transitions centred at 1.5, 10, 
16 and 49meV. In Ref. [20j, a weak feature near 70meV 
was also attributed to a magnetic transition, and the 
16 meV signal was found to be split into two peaks sep¬ 
arated by 2.5 meV consistent with earlier neutron^ and 
Raman22 spectra. The 10-16 meV transitions are clearly 
visible in Fig. 1 although they not resolved in this partic¬ 
ular data set<2£ The 49 meV transition is also present and 


Figure 2. (Color online) Neutron spectrum S(<j>,E) of poly¬ 
crystalline Tb 2 Ti 2 C >7 for fixed scattering angle <f>. The corre¬ 
sponding Q is given by h 2 Q 2 /2m n = Ei + E{ — 2(EiE{) 2 cos cj>, 
where E{ — Ei — E and m n is the neutron rest mass, (a) Spec¬ 
trum at three different average scattering angles recorded with 
Ei = 150 meV at T = 7K, showing magnetic peaks at 49 and 
61 meV and phonon peaks at 72, 91 and 99meV. (b) Com¬ 
parison of the low angle ({(/)} = 10°) spectrum at T = 7 and 
90 K, recorded with Ei = 65meV. 


confirmed as magnetic by the characteristic reduction of 
its intensity with Q due to the magnetic form factor. 
Above 50 meV there is a weak feature near 61 meV which 
decreases with Q, and there are peaks near 72 and 90- 
100 meV which increase with Q consistent with scattering 
from phonons. 
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The part of the spectrum with E > 20 meV is rep¬ 
resented in more quantitative detail in Fig. [21 Fig¬ 
ure (21a) shows the intensity as a function of energy for 
three different average scattering angles </> taken from the 
Ei = 150 meV run at 7K. The large peak at 49meV and 
the small peak at 61 meV both decrease with increasing 
(f> (i.e. increasing Q ), whereas the peaks at 72, 91 and 
99meV all increase with angle. This confirms that the 
latter three peaks are caused by scattering from phonons, 
to within the sensitivity of the measurement, whereas the 
49 and 61 meV peaks are from magnetic transitions. 

Figure (21b) compares <f> = 10° spectra at 7 and 90 K 
measured with Ei = 65meV. Two features stand out: 

(1) The 49meV peak in the 7K spectrum is broader 
than the resolution, with a width (full width at half max¬ 
imum) of 4.9 meV as compared with the instrumental res¬ 
olution of about 1.5 meV at this energy. The peak is also 
asymmetric, with a low energy tail that is not wholly 
accounted for by the asymmetry of the resolution func¬ 
tion. On warming to 90 K, the peak increases in intensity 
and shifts down in energy to 47.5 meV. Transitions from 
the thermally populated 1.5 meV first excited crystal field 
level account for some of the peak width at 7K but not 
all, and there are no sharp phonon peaks at this energy 
(see Fig. [TJ). These observations suggest that there exists 
an additional crystal field level at about 48 meV which is 
strongly connected to the 1.5 meV level. The peak posi¬ 
tions and widths at 7 K and 90 K would then be explained 
by transitions from the ground state and 1.5 meV level to 
levels at 48 and 49meV. 

(2) In the range 20-45 meV there is considerable addi¬ 
tional scattering at 90 K relative to 7 K. The 90 K scatter¬ 
ing does not display any particularly prominent features, 
but a small peak has grown at about 24 meV and there is 
a weak shoulder near 29 meV which coincides with a min¬ 
imum in the 7K spectrum. These temperature-induced 
features are not accounted for by the thermal population 
of phonons (the temperature factor for phonon scatter¬ 
ing [1 — exp(— E/k^T)]~ l increases by only 4% at 24meV 
on warming from 7 to 90 K). The additional scattering 
must therefore derive from transitions out of thermally 
excited magnetic levels, and since the separation between 
the peaks near 16 and 49 meV is significantly greater than 
25meV there must exist a hitherto undetected level be¬ 
tween these two. Indeed, thermally excited transitions 
from the 10 and 16meV levels to a level near 39meV 
would give rise to enhanced intensity around 24 and 
29 meV, as observed. 

In short, we have observed a magnetic peak at 61 meV 
which corresponds to a crystal field level not found in 
previous studies, and we have indirect evidence for addi¬ 
tional levels near 39 meV and 48 meV. 

To constrain the crystal held model for TTO as tightly 
as possible we shall also consider the transition intensi¬ 
ties. We determined the integrated intensities of the mag¬ 
netic peaks in the low angle spectrum at 7 K by fitting an 
asymmetric pseudo-Voigt line shape to the peaks. The 
parameters of the pseudo-Voigt function were determined 


Table I. Observed and calculated crystal-field transition ener¬ 
gies and integrated intensities from the spectrum of Tb 2 Ti 2 07 
measured at 7K. Numbers in parentheses are experimental er¬ 
rors in the last digit. Observed intensities have been corrected 
for the dipole form factor of Tb 3+ so that both I 0 b s and 7 ca i c 
are for zero Q. The 7 ca i c values include transitions from both 
0.0 and 1.4 meV levels with their respective thermal popu¬ 
lations at 7K. The best-fit crystal-field parameters used for 
the calculations are: Bq = 55.3, Bq = 370.4, 7?| = 128.0, 
Bg = 114.3, B$ = -114.3, Bt = 120.3 meV. 


Level 

Eobs 

(meV) 

E C alc 

(meV) 

7obs(Q=0) 

(mb sr -1 f.u. -1 ) 

7calc(Q = 0) 
(mb sr -1 f.u. -1 ) 

1 3 

0.0 

0.0 

- 

2918 

r+ 

1 3 

1.4(4£ 

1.5 

- 

2744 

rj 

10.2(4) 

10.3 

1740(60) 

1836 

1 't 

16.0(5^. 

16.1 

990(100) 

1010 

r+ 

1 3 

39(2) 

39.0 

1 

31 

r 2 

48(1) 

48.2 

230(35) 

62 

r t 

49(1) 

48.8 

J 

184 

r+ 

1 3 

61(1) 

60.8 

44(6) 

52 

r+ 

1 1 

- 

71.0 

- 

8 


a The first excited state is dispersive at low temperature with a 
band width approaching 1 meV ^9^14^20,25 
b Comprises two peaks separated by 2.5 meV 1/1,20 


according to an empirical implementation of the analytic 
line shape which contains a contribution from the veloc¬ 
ity selection via the chopper^ and an additional compo¬ 
nent due to the pulse width. The nonmagnetic (phonon) 
background was estimated from the high angle part of 
the spectrum, which was scaled to match the high en¬ 
ergy transfer part of the low angle spectrum. A list of 
observed energy levels and transition intensities at 7 K is 
given in Table [H For ease of comparison, the intensities 
have been extrapolated to zero Q via the Q dependence 
of the magnetic dipole form factor of Tb 3+ . 

The crystal field at the Tb 3+ site in TTO has point 
symmetry 3 m (D^d.) and is described by the Hamiltonian 

Hcf =BqCq + B*C£ + BUCt, - Ci) + B$C$ 

+ Bt(C% - Cl) + Bl(C e _ 6 + CD, ( 1 ) 

where B* are the crystal field parameters and are 
Wybourne tensor operators^ Diagonalisation of Hcf 
was performed in the intermediate coupling scheme with 
the program SPECTRE^ To speed up the calculation, 
the complete basis of the / 8 configuration of Tb 3+ (3003 
states) was truncated to the lowest 110 states (the com¬ 
plete 7 F term, plus the states 5 774, 5 Ge, 5 Lio and 
5 G 5 ) extending to 3.5 eV above the ground state. The di¬ 
agonalisation of Hcf included J-mixing within the trun¬ 
cated basis. 

The neutron spectrum of single-ion magnetic transi- 
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tions is given by2 

* i 

x|<r i |M_ L (Q) |r i )| 2 <5(£’ J - -Ei - £), (2) 

where (7ro/2) 2 = 72.7mb. The first summation is over 
the initial states T^ with thermal population pi, and the 
second summation is over the final states Tj. The Debye- 
Waller factor e~ 2W is taken to be unity at the low tem¬ 
peratures of the measurements. We assume the dipole 
approximation, in which case Mj_(Q) can be replaced by 
— f(Q)gjJ± where f(Q) is the dipole form factor, gj is 
the Lande g- factor and Jj_ is the component of the to¬ 
tal angular momentum perpendicular to Q. Calculated 
intensities are powder-averaged for comparison with the 
data. 

The six crystal field parameters in Eq. (©) were refined 
by a weighted least-squares fitting algorithm against the 
experimental data given in Table HI with the intensities 
expressed relative to the 10.2 meV peak. The squares 
of the experimental uncertainties were used as recipro¬ 
cal weights. The crystal field parameters determined for 
Ho 2 Ti 2 07 (Ref. |H) were used as starting parameters for 
the fit. The procedure converged to an excellent fit with 
X 2 = 0.95, where % 2 is the standard normalised goodness- 
of-fit statistic. Fits were also performed with different 
sets of starting parameters but no other acceptable dis¬ 
tinct solutions were found. In particular, fits starting 
from the parameters found by Zhang et did not con¬ 
verge. 

The best-fit parameters are given in Table U together 
with the calculated energy levels and intensities at 7K. 
The calculated values agree very well with the obser¬ 
vations, including the absolute intensities which have 
a systematic error of 5-10% from uncertainties in the 
vanadium calibration and in the corrections for attenua¬ 
tion and the Q dependence of the magnetic form factor. 
Figure [31 compares the predictions of the best-fit model 
with the spectra measured at 7 K and 90 K. Pseudo-Voigt 
line shapes have been used to model the resolution func¬ 
tion. The intensities are calculated in absolute units from 
Eq. ([2|) and have not been scaled to fit the data. Over¬ 
all, the agreement is very good. The main discrepancy 
is that the measured spectrum at 90 K is less structured 
above 20meV than the calculated spectrum, which sug¬ 
gests that the crystal field levels in this energy range 
broaden significantly with temperature. 

IV. DISCUSSION 

We have presented here a single-ion model based solely 
on measurements on Tb 2 Ti 2 07 that successfully de¬ 
scribes the general features of the observed magnetic 
spectrum of Tb 3+ in Tb 2 Ti 2 07 . The crystal field pa¬ 
rameters of our model are compared with previously pub¬ 
lished sets [converted to the Wybourne tensor parameters 


Table II. Comparison of crystal field parameters Bq (in meV) 
for Tb 2 Ti 2 07 from different analyses. The parameters are 
defined for the Hamiltonian in Eq. ©■ The effective Hamil¬ 
tonian in terms of Stevens operators O)?,— which applies for 
a basis comprising a single J level in the LS coupling scheme, 
is Rcf = JO 8kD k O k . The 8 k are reduced matrix elements 
which for Tb 3+ (J = 6) are 82 = —1/99, 84 = 2/16335 and 
8q = —1/891891. The parameters D k and Bq are related 
by Dl = \ q k Bq, where X° 2 = 1/2, A4 = 1/8, A| = y/35/2, 
Ag = 1/16, Ag = VT05/8 and Ag = v/531/16. 


Ref. 

to 

OtO 

So 

S 3 4 

So® 

Bl 

B% 

4 

53.6 

318 

146 

149 

-143 

67.6 

14 

60.9 

291 

103 

96.6 

-59.9 

97.5 

18 a 

56.0 

329 

95 

107 

-77.4 

109 

19 

67.3 

320 

119 

113 

-90.5 

101 

20 

144 

268 

162 

171 

349 

799 

29 b 

75.3 

329 

100 

111 

-79.6 

130 

This work 

55.3 

370 

128 

114 

-114 

120 


a The model in Ref. 11 d is a refinement of the models presented in 

Reis. Ql 03 

b Parameters for Ho 2 Ti 2 07 scaled with the point-charge 

relational Bq (Tb)/Bq (Ho) = (r fc )Tb/(r fc )Ho, where ( r k ) is the 
fc th radial moment of the 4/ electron distribution. 

of Hamiltonian (|T|) ] in Table HT1 Our parameters have sim¬ 
ilar magnitudes and signs to those obtained for TTO in 
Refs. @ an dEMl (which however relied on spectroscopic 
data on other rare earth titanates). The crystal field pa¬ 
rameters are also similar to those of Ho 2 Ti 2 07 (Ref. i29f ) 
and Pr 2 Sn 2 07 (Ref. H3) after taking into account the 
differences in the radial moments of the respective 4/ or¬ 
bitals. The lack of significant variation shown in the crys¬ 
tal field potential across different systems implies that 
the local structure and bonding is similar for different 
pyrochlore oxides. 

As deduced many years ago^i the low energy part of 
the spectrum at low temperature is composed of two non- 
Kramers doublets separated by approximately 1.5 meV. 
From our model, the largest components of the ground 
state and first excited doublet wave functions are found 
to be 

Fg"(0meV) = 0.968| 7 E 6 , ±4) =f 0.089| 7 F 6 , ±1) 

+ 0.110| 7 E 6 , t2) ± 0.181| 7 F 6 , t5) 
-0.083| 7 E 4 ,±4) 

r+(1.5meV) = 0.952| 7 F 6 ,±5) t0.176| 7 E 6 ,±2) 

+ 0.061| 7 E 6 , =Fl) ± 0.1941 7 W 6 , ^4) 
±0.134| 7 F 5 ,±5). 

The upper and lower signs give the two components of the 
doublet, and \ 2S+1 Lj,mj) identifies the spectroscopic 
term and rrij value for each component. The dominant 
components in the wave functions are from the Hund’s 
rule ground state term | 7 F 6 ), as expected, but there is 
a non-negligible admixture of the | 7 W 5 ) and | 7 W 4 ) states 
which in the free ion lie 250 and 410 meV, respectively, 
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Figure 3. (Color online) Comparison of the observed and calculated spectrum of polycrystalline Tb2Ti2C>7. The heavy solid 
(red) lines are calculated from Eq. ([2| with eigenstates derived from the crystal field parameters given in Tabic [U and the 
individual transitions are shown as thin (blue) lines. The broken (blue) lines are the estimated non-magnetic background, and 
the fitted elastic peaks centred on E = 0 are shown as pale (green) lines. 


above the ground state. For comparison, with the same 
crystal field model but working in the pure |~Jf^) basis 
with the Stevens operator form of T~Lcf instead of Eq. m 
we find that the overall splitting increases from 71 meV 
to 81 meV. A list of the wave functions for all the levels 
within the \ 7 Fq) manifold is given in the Appendix. 

Our results are consistent with most previous stud¬ 
ies, but are in contrast to the work of Zhang et al£l 
who found essentially the same two lowest doublets as 
above but in the reverse order, i.e. with \ 7 Fq, ±5) as the 
dominant components of the ground state doublet. This 
difference arises because Zhang et al. assumed in their 
analysis that the splitting of the 16meV peak was due 
to two distinct crystal field transitions. Although the 
model of Zhang et al. fits their own data well, our mea¬ 
surements reveal a number of problems with it. Firstly, 
the magnetic peak we observe at 61 meV was not iden¬ 
tified as a crystal field transition by Zhang et al., and 
therefore not used to constrain the model. Second, their 
model predicts a level at 101 meV with a very large cross 
section for transitions to it from the ground state dou¬ 
blet: 440mbsr -1 f.u. -1 at zero Q. This prediction is 
inconsistent with our data since such an intense transi¬ 
tion would produce a peak centred at 101 meV in Figs. |T| 
and ETa) which, after correction for the Tb ,1+ form fac¬ 
tor, would be of roughly the same size as the peak at 
49meV. Third, the model of Zhang et al. predicts a 


strong transition from the first excited doublet to a level 
near 72 meV. The zero-Q cross section for this transition, 
which should appear as a thermally excited peak near 
70meV, is about 270 mb sr -1 f.u. -1 at T = 90 K, incon¬ 
sistent with the 90 K spectrum shown in the lower right 
panel of Fig. [3] Finally, the significant intensity from 
thermally excited transitions observed between 20 and 
40meV — see Fig. Efb) — is not reproduced. 

The success of our model, which assumes that the split 
peak at 16 meV corresponds to a single crystal field level, 
raises questions about the nature of this excitation. One 
possibility is that the splitting could be the result of dis¬ 
order, e.g. Tb/Ti site mixing, which could produce two 
slightly different local environments for the Tb site. How¬ 
ever, the fact that no splitting is detectable in any of the 
other transitions, especially the F 3 doublets whose degen¬ 
eracy is lifted once the 3m symmetry is broken, suggests 
that there is only one Tb 3+ environment. Further, the 
splitting is too large to be explained by the same two-ion 
magnetic coupling model which accounts for the disper¬ 
sion of the first excited doublet.— Klekovkina and Malkin 
have suggested that the peak splitting is caused by cou¬ 
pling to a phonon with an energy near 16meV at the 
zone centre.- 18 Indeed, hybridization of acoustic phonons 
with the 1.5 meV and lOmeV crystal field levels has re¬ 
cently been observed in neutron spectra of TTOj2ii£ and 
strong thermally-induced phonon anomalies are observed 
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Table III. Eigenfunctions of the best-fit single-ion crystal field model for Tb 3+ in Tb 2 Ti 2 C> 7 . The crystal field Hamiltonian is 
given in Eq. 0, and the values of the best-fit crystal field parameters are given in the caption of Tabic [T] The list below is of 
the energy levels within the | 7 F§) manifold and gives only the coefficients of | mj) belonging to \ 7 Fq). A blank entry means a 
zero coefficient. 


mj 

0.0 

0.0 

1.5 

1.5 

Energy (meV) 

10.3 16.1 39.0 

39.0 

48.2 

48.8 

60.8 

60.8 

71.0 

-6 





0.145 

-0.179 


0.687 

0.678 



-0.033 

-5 

0.181 


0.952 



-0.182 




0.033 



-4 


0.968 


0.194 



-0.071 




-0.096 


-3 





0.686 

-0.675 


-0.134 

-0.173 



-0.066 

-2 

0.110 


0.176 



0.956 




0.057 



-1 


0.089 


0.061 



-0.058 




0.974 


0 






-0.105 



0.019 



0.982 

+ 1 

-0.089 


0.061 



0.058 




0.974 



+2 


0.110 


-0.176 



0.956 




0.057 


+3 





0.686 

0.675 


-0.134 

0.173 



0.066 

+4 

0.968 


-0.194 



-0.071 




0.096 



+5 


-0.181 


0.952 



0.182 




-0.033 


+6 





-0.145 

-0.179 


-0.687 

0.678 



-0.033 


and attributed to magneto-phonon coupling. 32 The exis¬ 
tence of magnetoelastic modes was suggested as an ex¬ 
planation for why frustration is not relieved in TTO by 
any conventional symmetry-breaking transitions at low 
temperatures!*^ Although our data cannot shed much 
light on this interesting question, we mention that the 
crystal field states in TTO do have strong quadrupole 
moments so a perturbation to the single-ion states involv¬ 
ing a coupling via orbitals is plausible. Measurements 
on single crystals are needed to identify any coupling to 
specific phonon modes and to see whether the size of the 
splitting varies throughout the Brillouin zone. 

As far as the physical properties of TTO are concerned, 
the most important result from this study is that the 
ground state is dominated by the components | 7 i*6i ±4). 
This means that at temperatures where the first excited 
doublet has negligible population TTO has strong Ising- 
like anisotropy with the moments confined to the quan¬ 
tization axis, i.e. to the local (111) directions. For the 
ground state doublet, we calculate that the components 
of the zero-field spectroscopic g-tensor parallel and per¬ 
pendicular to the quantization axis are g\\ = 10.7 and 
g l = 0 (g_i = 0.07 in a field of 1 Tesla), consistent with 
previous estimates^ 

V. CONCLUSION 

We have successfully refined a model for the crystal 
field in Tb 2 Ti 2 0y against the magnetic spectrum mea¬ 


sured by neutron spectroscopy. The model is in very 
good quantitative agreement with the experimental spec¬ 
trum, including the absolute intensity, and confirms that 
the single-ion magnetic properties at low temperature 
are controlled by two non-Kramers doublets separated 
by about 1.5 meV. The wave functions of these states are 
approximately | mj) = |±4) and |±5), respectively. Our 
analysis is more tightly constrained by experiment than 
has hitherto been possible, and resolves an uncertainty in 
the literature about the composition of the ground state. 
The splitting of the 16 meV peak remains unexplained, 
and needs further investigation. 
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Appendix A: Single-ion eigenfunctions of Tb 2 Ti 207 

The eigenfunctions for the best-fit model are given in 
Table urn 
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